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Single transient laser-induced microbubbles have been used in microfluidic chips for fast actuation of the 
liquid (pumping and mixing), to interact with biological materials (selective cell destruction, membrane 
permeabilization and rheology) and more recenty for medical diagnosis. However, the expected heating 
following the collapse of a microbubble (maximum radius ~ 10-35 /mi) has not been measured due to 
insufficient temporal resolution. Here, we extend the limits of non-invasive fluorescence thermometry using 
high speed video recording at up to 90,000 frames per second to measure the evolution of the spatial 
temperature profile imaged with a fluorescence microscope. We found that the temperature rises are 
moderate (< 12.8°C), localized (< 15 /mi) and short lived (< 1.3 ms). However, there are significant 
differences between experiments done in a microfluidic gap and a container unbounded at the top, which are 
explained by jetting and bubble migration. The results allow to safe-guard some of the current applications 
involving laser pulses and photothermal bubbles interacting with biological material in different liquid 
environments. 



Short lived vapor microbubbles can be created by the absorption of laser pulses in liquid environments. The 
rapidly expanding bubbles later collapse driven by the pressure of the surrounding liquid. The oscillation of 
the bubbles induce fast flow, intense shear, shock wave emissions and adiabatic heating. These physical 
mechanisms are utilized in many experiments with microfluidic chips and biological material, examples are: 
pumping^ mixing (fast vortices)^'^, dropplet coalescence^, single cell destruction^'^, cell membrane poration^'^ and 
microrheology^. Several of these applications rely on nanosecond laser pulses absorbed with the help of dyes^°, 
nanoparticles^^"^^, microparticles^'^^ and nonlinear absorption^'^'^^. Other biomedical laser applications include 
eye surgery^^, tattoo removaP^ and laser microdissection^^. Recently, photothermal nanobubbles were used for 
needle-free malaria detection^\ 

Despite the success of laser-induced bubbles in combination with microfluidic chips and biomedical applica- 
tions, there remains an open question about the amount and spatial distribution of heat generated by these 
microscopic bubbles inside confined liquid geometries. Experiments with macroscopic bubbles (radius ~ 
10 mm) created by the arrest tube method inside aim long tube have shown temperature increases of up to 
4°C during the collapse of the bubble^^. Experiments with microscopic laser-induced bubbles in combination with 
biological material hint at moderate temperature rises, since it has been possible to reversibly permeabilize and 
stretch biological cells while retaining viability^'^. Inside microfluidic channels, non- spherical bubble collapse may 
lead to liquid jetting^ which could result in enhanced cooling through thermal advection from the vorticity 
generated. 

To our knowledge there have been no measurements of the temperature distribution following the collapse of a 
laser-induced cavitation microbubble. Here we use fluorescence based thermometry^^"^^ combined with high 
speed recording to spatially and temporally resolve the temperature distribution after the oscillation of a single 
laser-induced microbubble inside a microfluidic liquid gap. The maximum radius (—10-35 ^m) of these bubbles 
is similar to bubbles used in some microfluidic experiments in combination with biological materiaF'^ and 
bubbles created with other photothermal methods In order to understand what are the effects of the geometry 
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Figure 1 | (A) Experimental setup. DM, dichroic mirror. HM, 50/50 half 
mirror. OBJ, microscope objective. Fl and F2, filters. (B) Geometry of the 
liquid containers, microfluidic gap and semi-infinite liquid. Also shows the 
plane where the temperature is measured and averaged over concentric 
circles with radius r centered where the laser is focused and where the 
temperature is the highest after the bubble collapse. (C) Temperature 
calibration: temperature increase AT as a fijnction of (Ird ~ 1). 

on heat deposition and dissipation, experiments are also done in an 
unbounded liquid (at the top) geometry to compare the results with 
those inside the microfluidic gap. 

Results 

Experimental setup. The experimental setup is depicted in Fig. lA. 
Single laser pulses (duration of 6 ns and wavelength of 532 nm) from 
a frequency doubled Nd:YAG laser (New Wave, Orion) are used to 
create the bubbles by focusing the pulses inside the containers 
(microfluidic gap and unbounded) filled with fluorophore solution. 
The beam is expanded and then coupled into an inverted microscope 
(Olympus IX-71) through a dichroic mirror (DM) that reflects light 
at 532 nm. Then the light is transmitted through a 50/50 half mirror 
(HM) and is focused at the bottom of the fluorophore container 
using a 40X (0.8 NA) water immersion microscope objective (OBJ). 

The experiments are done using two different containers for the 
fluorophore solution (Fig. IB): a thin microfluidic gap bounded at 
the top and at the bottom with #1 coverglass (0.13-0.16 mm thick) 
separated by spacers with a height of 19 ± 2 ^m; and a small petri 
dish with a #1 coverslip glass bottom filled up to a height of 3 mm 
(resembling a semi-infinite liquid), i.e. much larger than the size of 
the bubble. 

For the temperature measurements we use single-color intensio- 
metric thermometry, since it only requires a simple epifluorescence 
light path that can be combined with the optical elements that couple 
the laser pulses to the microscope objective (Fig. lA). Furthermore, 
the choice of single color fluorescence is justified under the assump- 
tion of a constant fluorophore concentration, which we will discuss 
later (Discussion). The fluorophore solution is excited by a Mercury 
fluorescent lamp (Fl: excitation filter for 480-550 nm) reflected by 
the HM into the microscope objective. The fluorescence is recorded 
by a high speed camera (SAl.l, Photron) at up to 90,000 fps for the 
semi-infinite liquid and 54,000 fps for the thin liquid gap. The fiuor- 
escence emission filter (F2: long pass 590 nm) is placed in front of the 
camera. 



Temperature calibration. In order to find the relation between 
fluorescence intensity 1 and temperature, a small glass container 
with a calibrated thermistor is filled with fiuorophore solution (see 
Methods) to a height ~3 mm and sealed at the top with a Peltier 
heating element. The container is placed on the same experimental 
setup and is slowly heated while the fluorescence intensity is recorded 
(100 frames for selected temperatures) with the high speed camera at 
a frame rate of 54,000 fps during several heating and cooling cycles 
(from room temperature at 21.6 to 65°C). The intensity for each 
frame is integrated and averaged over the set of frames. Photo - 
bleaching is not observed for times greatly exceeding the duration 
of the experiments (<7 ms). Additional data is shown in the 
Methods section. We find a linear relation between temperature 
and intensity: T(J) = a + hi in the explored temperature range. 
Hence, the temperature increase AT = T — 21.6°C can be 
expressed as AT = To(I^e/ — 1), where Tq is a constant and Irei=I/ 
J(21.6°C). Figure IC shows the results for six calibration experiments 
where the temperature increase AT is plotted as a function of Ird ~ 1- 
The linear fit to the data yields Tq = — 1 14.5 ± 1.6°C where the error 
represents the confidence bounds within 95%. 

Experiment. Single laser pulses are focused at the bottom of the 
containers (for both geometries), above the bottom coverslip 
(2 ^m). The molecules in the focal volume could be bleached by 
the focused laser light. An estimate for that volume can be 
obtained by considering the size of the ellipsoid at the focal 
volume using the diffraction limited spot size Wq = 0.41 ^m and a 
Rayleigh length of = 2. 16 ^m. This is equivalent to a volume of 0.7 
femtoliters, which is negligible compared with the volume of the 
liquid contained in the vicinity of the irradiated spot. These 
molecules should be mixed with the rest of the solution during the 
bubble dynamics. Hence if the molecules in the focal volume are 
bleached, we can neglect them. We also assume that the fluoro- 
phore concentration does not change after the cavitation event and 
the possible effects of the bubble dynamics and thermophoresis are 
negligible (Discussion). The energies used for the laser pulse are in 
the range of 1-9 ^J, resulting in bubbles with a maximum transverse 
radius {xy plane) between 10-35 ^m for both geometries. For each 
geometry the experiment is repeated five times at each laser energy 
and the temperature distributions T(r, t) are extracted for each 
experiment separately. 

Image analysis. To analyze the dynamics of the temperature 
distribution, we averaged several hundred images (300 for the thin 
gap and 600 for the semi-infinite liquid) prior to the laser irradiation 
which gives us the fluorescence intensity at room temperature. Then 
each of the next 300-600 frames recorded during and after the arrival 
of the laser pulse is divided by the averaged room temperature image 
(to obtain Irei). In this way we are left with 300-600 images of the 
relative fluorescence intensity (minus one or two saturated images 
that show the arrival of the laser pulse and the bubble dynamics). 
Figure 2 A shows representative Ird images (converted to AT) 
recorded at 54,000 fps before and after the arrival of a single laser 
pulse (8.6 ^J) in the microfluidic liquid gap. The images were filtered 
with a circular averaging filter with a radius of four pixels in order to 
increase the signal-to-noise ratio (Video file provided in 
Supplementary Video 1). The Irdi^^ t) images are analyzed 
without filtering and the relative fluorescence intensity is averaged 
over concentric circles (radius r on xy plane) around the irradiated 
spot in a range L = 100 ^m (Fig. IB). The number of points in the 
perimeter increases linearly with the radius, as a result, the standard 
error of the mean (1.96(7/v^, <J is the standard deviation and n the 
number of points) decreases with increasing r. Then, Irdir, t) is 
converted to AT(r, t) using the calibration. An example of AT(r, t) 
extracted from unfiltered intensity ratios for some of the cases in 
Fig. 2 A are shown in Fig. 2B. Reliable measurements of the 



SCIENTIFIC REPORTS | 4 : 5445 | DOI: 1 0.1 038/srep05445 



2 



-28.1 las ^^g^H 27.4 |Lis 64.5 jus 138.5|as 

231.1|as 323.7 |Lis 416.3 |as 508.9 |as 601.5 |as 

m m ♦ 

694.1 iLis 786.7 iLis 879.3 |lis 971.9 |lis 1.064 ms 




27.4 las 



231.1 las - 



508.9 las 



0 20 40 60 80 100 

r(|Lim) 

Figure 2 | (A) Temperature distribution filtered with a circular average to 
reduce noise (thin gap, Eiaser = 8.6 fij). Video file provided in 
Supplementary Video 1. Size of each image is 101 X 116 jum (W X H). 
(B) Temperature distribution averaged at every r with no filtering. The 
dotted line represents a fit to a simple heat diffusion model. 



temperature distributions are made at —35 fis after the arrival of the 
laser pulse, that is, the first measurement is taken at 25 to 33 fis after 
the collapse of the bubble. 

Figure 3 compares the maximum temperature rise, distribution, 
bubbles dynamics, and cooling rates for the two geometries (black 
circles: microfluidic gap, red triangles: semi-infmite liquid). First, the 
maximum temperature rise above the initial temperature distri- 
bution ATo(r = 0, to) is plotted in Fig. 3 A as a function of the laser 
energy. The maximum temperature is higher for the microfluidic gap 
for similar laser energies. The largest measured temperature rise is 
11.9 ± 0.9°C and 8.9 ± 1.3°C for the microfluidic gap and semi- 
infmite liquid respectively. If we consider that a similar amount of 
heat is deposited in both geometries, less liquid in the vicinity of the 
irradiated area for the microfluidic gap could explain the higher 
temperatures. The measured radial width (half width at half max- 
imum) of the initial temperature distribution as a function of the 
laser pulse energy is plotted in Fig. 3B. The initial widths are similar 
for both containers. Figure 3C shows the lateral {xy plane) maximum 
bubble radius which was captured in another set of experiments with 
similar energy range of the laser pulses but recorded at 450,000 fps 
and imaged in bright field. A bubble created inside the microfluidic 
gap is shown in Fig. 3C (inset). The maximum radius squared has a 
linear relation to the laser energy R^^^ = aEiaser and it is very similar 
for both geometries. The constant a for the thin gap a2d — 152.8 ± 
18.4 lim^lpi], is slightly larger than for the semi-infmite hquid con- 
figuration Uo^d = 133.7 ± 5.8 fim^/fi]. 

Discussion 

In the semi-infinite liquid container the bubble has a hemispherical 
shape, hence we can use the 3d Rayleigh-Plesset equation to calculate 

the oscillation time Tqsc = ^-^^RmaxV p/ (poo —pv)^ where p is the 
liquid density, poo is the hydrostatic pressure (10^ Pa), and py is the 
vapor pressure (2330 Pa for water). For the largest hemispherical 
bubbles the oscillation time t^sc is under 7 ^s, consistent with the 
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Figure 3 | (A) Maximum temperature rise A Tas a function of the laser pulse energy (E/^^g^) for both chamber geometries. (B) Halfwidth at half maximum 
of the initial temperature distribution after the bubble collapse (Wq). (C) Maximum bubble radius R^ax as a function of laser energy obtained from 
experiments in bright field at 450,000 fps. The dashed lines, ^/ aEiaser > where a is the constant. Inset: sample of high speed video frames of a single cavitation 
event (microfluidic gap). (D) Time it takes the maximum temperature to drop by 1/e (t) as a function of Eiaser- Thin liquid gap: black circles, semi-infinite 
liquid: red triangles. 
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bright field high speed recordings. The bubbles created in the micro- 
fluidic gap are not planar, since R^ax is on the order of the height H of 
the gap. For these bubbles Tosc increases as the ratio ?] = H/R^ax 
decreases^^. High speed recordings in bright field showed that Tosc 

< 10 ^s. In the limit ^ ^ 0, the shrinkage of the bubble is planar and 
can be fitted to a 2d Rayleigh-Plesset equation. There are different 
collapse scenarios that depend on rj^'^. For rj in the range of (0, 0.4) 
there is planar collapse between both boundaries, upward jetting 
towards the upper boundary is expected for (0, 1), while for rj within 
(1, 1.4) the bubble collapses neutrally. Finally, for rj > lA jetting is 
directed toward the bottom boundary. In the microfluidic gap R^ax 
are between 15-32 ^m, hence rj has a range between 0.5 and 1.5, 
which should result mainly in upward jetting for the highest R^ax or 
Eiaser ^nd neutral collapse and probably downward jetting for the 
smallest bubbles. 

The data analysis on the dynamics of the temperature distribu- 
tions are done under the assumption that the fluorophore concen- 
tration does not significantly change after the cavitation event, that 
after the collapse of the bubble there are no appreciable gradients in 
the fluorophore concentration set by the bubble dynamics or by 
thermophoresis. First, we consider the diffusion coefficient of the 
fluorophore molecules (10 kDa TMR dextran at 10 mg/ml) in water. 
The diffusion coefficient as a function of molecular weight for FITC 
(fluorescein isothiocyanate) dextran solutions between 4 and 
2000 kDa can be calculated using D = oc (molecular weight)^ (cmV 
s), a = 2.71 X 10"' for dilute solutions (0.1 mg/ml) at 25°C°. The 
calculated diffusion coefficient for our 10 kDa sample is 8.97 X 
10"^ cm^/s ~ 90 fim^/s. The effect of the high concentration 
(10 mg/ml) in our TMR solution is to decrease the diffusion coef- 
ficient by about 7%^\ However, we can use the higher value of D = 
90 fim^/s as an upper bound, which is much smaller than the thermal 
diffusivity of water (1.5 X 10^ fim^/s). Next, we consider the effects of 
the bubble dynamics. During expansion it is possible that there could 
be an increase in fluorophore concentration around the edges of the 
bubble. However, in the collapse there is jetting and vorticity that 
should quickly mix the solution. The longest lifetime of the bubbles is 
Tosc ~ 10 ^s and our earliest measurements are done at 25 to 33 ^s 
after the collapse of the bubble. The fluorophore concentration after 
the bubble collapses should be very similar to that before the bubble 
generation. Another possible source for a concentration gradient is 
the thermal gradient during absorption of the laser energy and in the 
bubble collapse, which could cause a redistribution of the fluoro- 
phore by thermophoresis. This effect is estimated using reported 
Soret coefficients (St) for aqueous solutions of dextran with similar 
concentrations (1, 5 and 10 mg/ml) and higher molecular weight 
(FITC dextran 86.7 kDa, D ~ 2.8 X 10"" cmVs = 28 firnVsY^ which 
underestimates the diffusivity in our sample by a factor of three. The 
reported St^ at the temperature range of interest for this experiment 
(21-35°C, after the bubble collapse) lies between -0.05 and -0.025. 
At higher temperatures (~ 45°C) Sr vanishes and becomes positive 
(migration toward colder regions) and at 55 °C, St ~ 0.025. In our 
experiment, the timescales for energy deposition and heating are very 
short: 6 ns during absorption and a submicrosecond bubble collapse. 
The initial heating during absorption of the laser pulse (6 ns) is 
followed by cooling during expansion of the bubble and finally heat- 
ing with the collapse of the bubble. The highest temperatures should 
be achived during absorption of the laser pulse. However the very 
brief timescales of the processes allow us to neglect thermophoresis 
even if St increased by several orders of magnitude due to high 
temperatures in our sample. The fluorophore diffusion length in 
10 ^s (upper bound) and considering St = 250 (4 orders of mag- 
nitude the value at 55°C) is 2^ (SrD)^^300 nm. Similar results are 
obtained if we assume that the initial thermal gradient (T ~ (21, 
35)°C and St ~ (-0.05, -0.025)) set after the bubble collapse is 
constant (continuous heating) for the duration of the experiment 

< 7 ms. Hence we can assume that the fluorophore concentration 



is constant which justifies the use of single color intensiometric ther- 
mometry and allows us to analyze our data under the assumption of a 
constant thermal diffusion coefficient or thermal diffusivity. 

Next, we look at the heat diffusion and re-thermalization of the 
liquid. The initial temperature distribution (AT(r, t) measured just 
above the bottom boundary) should decay exponentially in time 
according to the heat diffusion equation aV^T = dfT, where a is 
the thermal diffusivity. To simplify the notation we refer to the 
measured AT(r, t) as T(r, t), which can be expressed as a product 
of a spatial and a temporal part T(r, t) = F(r)G(t). The temporal 
solution is a decaying exponential exp( — lt), where 1 is proportional 
to the diffusivity a. To get an idea of the decay rate 1 we measure the 
time that it takes for the initial maximum temperature T^axi^ = 0, ^o) 
to drop to a value of T^^Je, we label this time by t. t is plotted on 
Fig. 3D as a function of the laser pulse energy. The mean value of all 
the data points for the thin liquid gap is T2d = 326 ± 26 ^s (central 
black line bounded by standard error of the mean) and for the semi- 
infinite liquid is = 412 ± 36 ^s (red line). Interestingly we 
observe that for the semi-infinite liquid all the data points are close 
to the shaded area within the error bars, suggesting that T^d (or 1^^) is 
constant or has little change for different Eiaser^Rmax- In the case of the 
microfluidic gap the variation of T2d hints at different dissipation 
rates that depend on Eiaser^Rmax which can be related to different 
collapse scenarios (rj) that can include jetting. The time for the tem- 
perature rise to drop to a 0.5% is about 3t which gives a time of 1.3 
and 1 millisecond for the semi-infinite liquid and microfluidic gap 
respectively, consistent with our measurements. 

To examine the effect of the boundaries we calculate the thermal 
diffusion length during the short times of the process (~ 1 milli- 
second). The diffusivities for water and glass are 1.5 X 10^ and 3.4 X 
10^ fim^/s respectively. For water the thermal diffusion length over a 
time of 1 ms is 24.4 fim, while for glass it is 36.8 fim, which are 
smaller than the thickness of the cover glass. 

Model. The evolution of the measured temperature profiles are fitted 
to the analytic solution of the heat diffusion equation (neglecting the 
boundaries) in order to extract effective thermal diffusivities. For the 
thin gap we treat the heat diffusion as a 2 dimensional (2d) problem 
because the transverse range {L= 100 ^m) in T(r, t) is larger than the 
height of the gap H. The case of the semi-infinite liquid is treated as 
an isotropic 3d temperature distribution (even though we only have 
half a sphere) problem in spherical coordinates. For both problems 
the boundary condition is T(r = L, t) = 0. The transverse range L is 
chosen because for the range of laser energies used, we do not detect 
temperature increases at that distance. Even if the fitted values of the 
extracted effective diffusivities changed by choosing a different L, we 
are interested in the relative differences between the extracted values 
for both containers. The Fourier series solution for the heat equation 
(with no sources) in 2d^^ and 3d^^ are: 
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The amplitudes are fitted to the initial temperature profile T(r, ^o) 
and the diffusivity is fitted to a temperature profile at a later time 
(Dotted lines in Fig. 2B). The predicted evolution of the temperature 
profiles are in good agreement with the measurements. The fitted 
effective diffusivities a for each experiment (single bubble) are shown 
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Figure 4 | (A) Thermal diffusivity (a) plotted as a function of Eiaser- Thin gap: black circles, semi-infinite liquid: red triangles. (B) Microfluidic gap: a as a 
funtion of rjy oc extracted from fits to 2d model (circles) and 3d model (squares). Drawing: expected collapse and jetting conditions for different ranges of rj. 



in Fig. 4A. Considering all the data the effective thermal diffusivities 
area = (2.51 ± 0.21) X 10' ^mVs and (1.34 ± 0.20) X 10' /imVsfor 
the thin liquid gap and the semi-infinite liquid respectively (black 
and red lines in Fig. 4A). In the semi-infinite liquid configuration we 
recover an effective thermal diffusivity very close to the value of water 
and all the data is centered around the mean value for all the pulse 
energies, similar to t^^ in Fig. 3D. This suggests that the temperature 
distribution is isotropic (in half a sphere). In contrast, for the micro- 
fluidic gap the fitted values of a are scattered around the mean value 
and depend on the laser energy (similar to T2d) hinting at a 3d 
temperature distribution with different dissipation rates that depend 
on Eiaserl^max ^^d hence the different collapse scenarios given by r\. 
Furthermore, the effective average value of the thermal diffusivity is 
almost twice that of water (from the ratio i2dl'^3d — ^idl^^d ~ 1-3 we 
expected an average value of a.2d ~1.3 X as^). The cylindrical distri- 
bution of the 2d model overestimates the average of the thermal 
diffusivity. 

Figure 4B shows the thermal diffusivity data for the microfluidic 
gap as a function of r\, where the black circles and red squares are the 
extracted as from the 2d and 3d model respectively. The average 
value of the 3d fits is a = 1.27 ± 0.15 X 10' ^mVs, which under- 
estimates the expected average value of d^d (1.3 X a^^. Both data sets 
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Figure 5 | Relative fluorescent intensity at room temperature recorded at 
54,000 fps. Room temperature data (bold line) recorded for about 0.37 s 
(upper time scale), each point corresponds to the average of one frame, 
adjacent points separated by 100 frames. Thin lines are samples of the 
calibration data with a duration of 1.9 ms (lower time scale). 



show essentially the same trend suggesting variable dissipation rates 
that depend on with the faster dissipations for the cases with 
upward jetting (0.5 < rj < 1), and the slowest for the conditions of 
neutral collapse and downward jetting which is closer to the case of 
the collapse on the semi-infinite liquid geometry {rj oo). 
Interestingly, we observe in both data sets that the value of a 
decreases for ?] approaching 0.5, which hints at a lower heat dissipa- 
tion rate for the case of a slower planar collapse (0 < < 0.4). Finally, 
we can estimate the amount of heat deposited after the collapse of a 
bubble by integrating the initial temperature distribution (assume 
isotropic on half sphere) multiplied by the specific heat, which yields 
between 5 and 10% of the laser pulse energy. 

In summary, we have shown that for the range of laser pulse 
energies and maximum bubble sizes, the temperature rise resulting 
from the oscillation of a single bubble is moderate, short lived and 
localized. Interestingly, the temperature distribution is a function of 
the geometry bounding the flow, which we explain with different 
jetting conditions (heat dissipation rates) during the collapse of the 
bubble. These results allow to safe-guard some of the current appli- 
cations involving focused laser pulses and photothermal bubbles. 

Methods 

Fluorophore solution. Water soluble tetramethylrhodamine attached to 10,000 MW 
dextran (TMR-Dextran) acts as the molecular thermometer. The fluorophore 
solution was prepared by dissolving the TMR-Dextran powder (Invitrogen, D1816) 
in distilled water (10 mg/ml). Even at the fastest recording times of 90,000 fps, the 
exposure time is much longer than the relaxation time of TMR which is on the order 
of nanoseconds. 

Figure 5 shows the relative fluorescent intensity recorded at 54,000 fps for more 
than 0.35 s (bold line, upper timescale), which greatly exceeds the duration of each 
experimental run. Hence we can discard photobleaching due to the excitation lamp. 
Also in Fig. 5 we show some of the data for the calibration (100 frames, lower 
timescale). 

Supplementary Video 1. Video for the case shown in Figure 2 (evolution of the 
temperature distribution), 102 X 102 /im (W X H) filtered with a circular averaging 
filter with a radius of four pixels. The total duration is 6.5 milliseconds. 
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